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Abstract 

In this paper we study parametrized description of neutrino mixing from a phenomenological 
point of view. Motivated by recent experimental results, we consider the case where a leading 
order mixing matrix receive small but non-negligible higher order contributions. A perturbation 
method to incorporate subleading contributions and its applications to tri-bimaximal mixing are 
discussed. We show that mixing matrices similar to tri-bimaximal-Cabibbo mixing can be obtained 
by straightforward choices of perturbation parameters. To achieve a better agreement with the ex- 
perimental data without increasing the number of free parameters, we impose a simple phenomeno- 
logical relation from which a trimaximal-like mixing matrix parametrized by U e 3 = s'm9i3e~ tip can 
be derived straightforwardly with no additional requirement on the perturbation parameters. It 
can describe the current global fit to three neutrino mixing to a good accuracy. Its theoretical 
explanation and phenomenological applications are discussed briefly. 
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I. INTRODUCTION 



Neutrino mixing is one of the most extensively studied topics in neutrino physics, which 
to a large extent is due to the tempting mixing pattern observed in neutrino oscillation 
experiments. It provides a clear evidence implying a non-trivial but perhaps simple flavor 
structure of the lepton sector. As a consequence, many interesting mixing patterns are 
proposed, including tri-bimaximal mixing (TBM) [1], bimaximal mixing [2|, golden-ratio 
mixings [3]], democratic mixing and hexagonal mixing j^j, etc. They often appear as 
exact mixing matrices, e.g. in models where lepton mixing is determined by a underlying 
discrete flavor symmetry (see, e.g. [6|). By exact mixing matrix we mean the mixing matrix 
that does not depend on other parameters including lepton masses 1 . For recent works and 
other exact mixing patterns, see e.g. (9]. 

This paradigm works quite well before the recent discovery of a relative large #13 (for 



recent global fits, see 



10Nl2|). which signals the deviation from those exact mixing patterns. 



For example, a recent global fit by M. C. Gonzalez-Garcia et al. [12] gives 



sin 2 fl 12 = 0.30 ±0.013, 



sin 2 9 



23 



0.41 



+0.037 
-0.025 



0.59; 



^0.021 
-0.0225 



(1) 



sin 2 # 13 = 0.023 ±0.0023. 



where two la ranges for sin 2 823 are given because the present data cannot resolve the 623 
octant degeneracy 12j . 0y are the mixing angles defined in the standard parametrization. 
Using the global fit values given in Eq. ([T]), one can evaluate the squared absolute values of 
the mixing matrix elements \{U u )ij\ 2 , which can be written collectively in a matrix as 



( 0.684^ 0.293i°;™ 



IK** 



v)i] I 



0.401 



V 



0.023 

f0.036 , 
-0.024 t 

+0.024 , 



^0.0023 
-0.0023 
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(2) 
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CP phase whose experimental value has a relatively large error. These omitted elements 
are often used to derive relations between mixing parameters including the CP phase, but 



0.036 

v)%3 \ with j ^ e and j 7^ 3 are omitted because they are affected by the Dirac 



The corresponding effective neutrino mass matrices are sometimes calied form diagonalizable matrices 

Bfl. 
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in our discussion they are determined by other parameters and hence the CP phase can be 
extracted by the standard formalism. 

Motivated by the experimental results, many variants of those exact mixing matrices 



are proposed (see, e.g. |l3Hl5l|). In many cases neutrino mixing can be described by the 
deviations from an exact mixing matrix. The later corresponds to the dominant leading 
order (LO) contributions determined by underlying flavor dynamics. The former is induced 
by subleading contributions including radiative corrections, charged lepton corrections, etc. 
(see e.g. [3]). In general, the deviations from a given mixing pattern can be described by 
the corrections to mixing angles and the Dirac CP phase 2 , which can be studied in concrete 
models or in model independent ways. 

In this paper, we discuss parametrized description of neutrino mixing from a phenomeno- 
logical point of view. We treat subleading contributions as perturbations and a perturbation 
method to incorporate them will be discussed in Sec. 2. In this method, perturbations act 
directly on the mixing matrix. The perturbed mixing matrix can be expressed in terms of 
linear combinations of the column vectors in the LO mixing matrix. In some cases it is 
simpler than dealing with mixing angles, and the physical relevance of the perturbations 
introduced in this way is more transparent. Detailed discussions will be presented in Sec. 2. 
Although this method can be applied to more general cases, we concentrate on its appli- 
cations to tri-bimaximal mixing. The derivation of two mixing matrices that are slightly 
different from tri-bimaximal-Cabibbo (TBC) mixing [l^] is discussed in Sec. 3. For relevant 



works, see e.g. 18|, Il9] and other references given in this paper. 



We notice that the original TBC mixing and the TBC-like mixing patterns agree with the 
data only marginally. To achieve a better agreement, sizable corrections to #12 and #23 are 
also needed and it seems necessary to introduce more parameters. In general, more param- 
eters mean more flexibility and more chances to survive the future data, but a description 
with less parameters can lead to simplified descriptions of neutrino phenomenologies and 
may provide more clues to underlying physics. In this respect, besides model dependent 
theoretical exploration, it is also worthwhile to look for and study simple descriptions sug- 
gested by the experimental data. Phenomenological or empirical relations can provide some 
good examples, e.g. the relation between 13 and the Cabibbo angle 6 C : 13 ps 6 c /V^ Q], 



2 Majorana phases can be ignored since they do not affect neutrino oscillations ■ 
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the quar k-lep ton complementarity 
relations 



«j|,-™P^Wy, bi-.axge mixing 



24j . bi-pair mixing relations |25j], etc. Besides their theoretical implications, they 
may also give rise to economical but rather accurate descriptions of neutrino mixing and 
phenomenology. 



In this paper, we introduce a relation, i.e. | (U v ) x | 2 — | (U, 



2/3, which agrees with 



the data quite well, as can be seen from Eq. ([2]). Detailed discussion about this relation is 
presented in Sec. 4 at where it is shown that this relation leads to a mixing matrix that 
can describe the current global fit to a good accuracy, including sin #23 . It is derived by 
the method discussed in Sec. 2 and parametrized by U e % = sin^e - * 1 ^. Hence it can be 
regarded as an improved TBC mixing. Moreover, it is also a trimaximal-like mixing matrix 
with the trimaximal condition |l5j being perturbed by a small correction. Hence we refer 
it as trimaximal- Cabibbo mixing. Its resemblance to trimaximal mixing suggests a possible 
explanation to itself and the relation from which it is derived, which will be discussed briefly 
in Sec. 5 where a brief summary and its phenomenological applications are also presented. 



II. PERTURBATIONS TO MIXING MATRIX 



Below we discuss a perturbation method to incorporate subleading contributions into a 
LO mixing matrix, which is useful when perturbations or corrections are small. We begin 
with writing the lepton mixing matrix U v as 

U U = U°T or U u = TUl 

where T is a unitary matrix describing the perturbations to U®, which is a LO mixing matrix. 
We concentrate on the case where the perturbation matrix T multiplies U® from the right, 
i.e. U u = U®T. The formalism developed for this case can also be used for the case where 
U u = TU®, which will be discussed briefly at the end of this section. 
U u and U® can be written as 

U U = (U 1: U 2 ,U 3 ), U® — (Ki,K 2 ,K 3 ). 

The unitary matrix T can be parameterized by three angles denoted by and six phases. 
As usual, angle corresponds to the rotation angle in the plane. First consider the 
simplest case with two vanishing angles, i.e. perturbation by a single rotation. For example, 
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in the case where s± 2 = £23 = and ^13 7^ 0, T is given by 



^ cos 9l 3 sin 9L 3 e %a ^ 



T 



\ 



1 

sin<^ 13 e m cos ^3 



where = diagje^ 1 , 1, e luJ2 }. Since P w and e luJ can be adsorbed by Majorana phases or 
charged leptons, we will ignore them in the following discussions. It then follows that 



where 



U x = (K, - x*K 3 )/f u U 2 = K 2 , U 3 = (K 3 + xK^/h 



fx = h = \A + \ x \ 2 ' x = tan ^i3e" 



(3) 



Since x is an arbitrary complex number, instead of angles and phases, below we will use 
parameters x and y (see below) as perturbation parameters because they can simplify the 
expressions when |x| and \y\ are small. 

The above procedure can be applied consecutively and iteratively. Below we consider two 
cases that are relevant to later discussions. In the first case £7° is perturbed by a rotation 
in (1,3) plane followed by a rotation in (1,2) plane. By applying Eq. twice on these 
rotations, it is straightforward to find that 



U 1 = (K 1 -y*f 3 K 2 -x*K 3 )/f 1 , 
U 2 = (f 3 K 2 + yK 1 -x*yK 3 )/f 2 , 
U 3 = (K 3 + xK 1 )/f 3 , 



(4) 



where 



A = / 2 = V(i + N 2 )(i + M 2 ): h = VT+W- 

For small \x\ and \y\, higher-order terms can be ignored and hence 

m-il-a) K 1 -y*K 2 -x*K 3 , 
U 2 ~(l-b)K 2 +yK 1 -x*yK 3 , 
U 3 ~ (l-c)K 3 + xK u 



(5) 



(6) 



where 



(|x| 2 + M 2 ) /2, b=\y\ 2 /2, c=\x\ 2 /2. 
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(7) 



As another example, we consider the second case where a (1, 3) rotation followed by a 
(2, 3) rotation. Similarly, one has 

U x ~{l-c)K x -x*K Zi 

U 2 ~{1- b) K 2 - y*K 3 - xy*K h (8) 
U 3 2i(l-a)K 3 +yK 2 + xK 1 , 

Other mixing matrices associated with two or even more rotations can be obtained in 
the same way. This method can also be used in the case where the perturbation matrix T 
multiplies U® from the left, i.e. U u = TU® by applying this method to its transpose, i.e. 

Uu = u° u f. 

Obviously, the perturbation method discussed above is different from the one dealing 
directly with mixing angles. It leads to a simple formalism for the construction of perturbed 
mixing matrices. It can also lead to some interesting findings, e.g. the trimaximal-Cabibbo 
mixing derived in Sec. 4. This method can be applied to any LO mixing matrix, no mat- 
ter how it is parameterized. Note that even if the LO mixing matrix is in the standard 
parametrization form, the resulting perturbed one may not. This is not physically signifi- 
cant, although it may require slightly more work to extract mixing parameters used in the 
global fits, especially the CP phase. In addition, in some cases, the physical relevance of 
the perturbations introduced in this way is more transparent. For example, in the case 
discussed just above the perturbation matrix can be related to charged lepton corrections. 
Trimaximal-Cabibbo mixing is also an example, which will be discussed later. 



III. TRI-BIMAXIMAL-CABIBBO MIXING 



In the rest of this paper we consider perturbations to tri-bimaximal mixing, i.e. 

( 2 y/2 \ 



t^TBM — 



1 

\7f5 



-1 V2 -V3 

It may be instructive to show some simple applications first. We begin with the simplest 
case discussed in the previous section. From Eq. ([3]), one finds that the second column 
remains intact. Hence the perturbed mixing matrix U u is a trimaximal mixing matrix (see, 



e.g. 



15j). Since this case is very simple, below we consider the other two cases. 



In the follo wing discussion we denote U e % by \e~ lLp / y2, as in tri-bimaximal-Cabibbo 

(TBC) mixing 14] which is also parametrized by sin #13 whose global-fit value is in good 

agreement with the relation sin 13 = Xc/y/2 where \ c = 0.2253 ± 0.0007 [26J is the Wolfen- 

stein parameter. We consider first the case where U v is given by Eqs. flH]) from which it 

follows that _ 

1 [2 A 

U rf • 'V 3 =T2 

where terms with order 0(X 3 ) or higher are neglected. Perturbation parameters x and y 
can be written as _ 

x = x \<* p ^, y = W)Ae^y|. (9) 

Then one has 

1 (10) 



\ya + xo\ = 1 

Because U v contains only one physical Dirac CP phase, for simplicity we require that xq 
and yo are real. In addition, we require that they do not depend on other parameters. In 
principal, one may use any x and yo as long as \x\ and \y\ are small. As an example, let 
yo = 1/3 and xo = 2/3. Then from Eqs. (jHJ) and (jHJ) one has 



U„ 



2 
3 



1-^ 

1 6 

- \e iv - 



y/3 \ l 12 



V2 



\ V6 



6 
12 



A i<f 



+ 



1 



1 - 



which is similar to the TBC mixing introduced in jl4|. 
U u are given by 



/ 



C(A 3 



V2 V ' 1 
The mixing angles extracted from 



sin 2 fl 



12 



- - '—, sin 2 6*23 



2' sm = y 



1 _ A_ 2 
3 9 

Substituting Ac for A, one finds that the deviation of sin 2 6*12 from its TBM or TBC value 
is A^/9 = 0.0056, which is the only difference between the mixing derived above and TBC 
mixing and much smaller than the experimental error. 

As another example, consider the case where U u is given by T[/tbm or its transpose 
U v = UtbmT. We use Eqs. (JH]) to construct U u . Denote £/tbm by (Ki, K 2 , K3) and substitute 
Ki into Eqs. ([6]), from |C/ 3e | = A/v2 one finds that \y* — x*\ = A. As above, one may let 
x = x Xe tip and y = yo^e lip and then it follows that \y — x \ = 1. Letting y = —x = 1/2 
leads to another TBC-like mixing matrix with mixing angles given by 



sin 2 # 



12 



-, sin 2 fl 



23 



A 2 

- — , sin 2 0i3 



A 2 
2 ' 
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IV. TRIMAXIMAL-CABIBBO MIXING 



As we can see, form the TBC-like mixings discussed above, one has sin 2 9 12 — sin 2 9j 2 BM = 
1/3 and sin 2 6» 23 ~ sin 2 9^ M = 1/2, which fit the data within the range between the 2a 
and 3a experimental bounds. To achieve a better agreement with the data, corrections to 
9 12 and #23 should also be taken into account. As discussed in the introduction, to avoid 
introducing more parameters, we will impose a phenomeno logical relation. We note that, 
to some extend, several well-known mixing patterns such as TBM and BM can be regarded 



as phenomenological mixing patterns. Trimaximal mixing [15] and TBC mixing [14] can 
also be derived with certain phenomenological relations in mind. The relation used in our 
discussion is given by 

l(^) el | 2 - l(^) e3 | 2 = § (11) 

which is well consistent with the present data. Like most other phenomenological or empir- 
ical relations, it can hardly be generated as an exact relation. Note that this relation is also 
satisfied by TBM and hence can be considered as a phenomenological constraint on the de- 
viations from TBM induced by higher-order corrections. Its possible theoretical explanation 
and phenomenological applications are discussed in the next section. 

We discuss first its implication on mixing parameters. As in the previous section, we 
denote |(t/j,) e3 | by \/\/2 and then from Eq. (TTTT) we have 

l(0iOeil a = f + y> l(^) e2 | 2 = ~-A 2 , l(tU 3 | 2 = y (12) 

from which it follows that 

sin 2 £ 13 = y, sin 2 # 12 ^(l-^A 2 

where higher order terms are neglected. Substituting the global-fit value of sin $13 given 
in Eq. (JT|) into the second equation, one finds that the deviation from sin 2 9j 2 BM = 1/3 is 
5 sin 2 #13/3 = 0.038 ± 0.0038, which is rather large when compared with the experimental la 
error of sin 2 9 12 ■ Therefore, to satisfy Eq. (JTTTl . the correction to 9j 2 BM cannot be neglected. 
Also note that this relation does not constrain the atmospheric mixing angle #23- Never- 
theless, below we will show that, even without making a particular choice of perturbation 
parameters, the perturbations described by Eqs. ([6]) lead straightforwardly to an improved 
TBC or trimaximal-like mixing matrix that can fit the data very well, including #23 . 
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We now discuss the perturbations to TBM satisfying Eq. (ITT]) . We begin with the case 
where £/tbm is perturbed from the right, i.e. U v = UtbmT. Since \{U v ) el \^ > 2/3, it must 
receive a contribution from the first element of Kj BM , the second column of U^bm- Therefore 
we use Eqs. to construct U u , which is the only choice we made. From the last equation 
in Eqs. §6§ one finds that \{U U )J = X/V2 leads to \x\ ~ \/3\/2 and, from the first equation 
and Eq. (fl2|) one has 

2 



1 



\x\ 2 + \y\ 2 



1 



y 



2 A 2 

3 + y 



For simplicity, y is assumed to be real. Since \x it is straightforward to solve the 

above equation for y. Ignoring higher order terms, one has 



x 



73 



Ae 



-up 



3V2 



A 



(13) 



Then from Eqs. ([6]) it follows that 



( 



\ 



3Ai 
2 



Ae~^ 
V2 



!(! + ¥) 71 

3Xe^ _ 15\ 2 \ 1 / -i I 3AM A< 

2 8 i V3 I 1 4 J v/2 1 1 _f ' 



X ( ~\ _i_ 3Ae' 



3A 2 



2 8 

15A 2 \ 1 A i 3\ 2 \ 1 ( -l _ Xe-^ _ 3A2_\ 



+ 0(X 3 ). (14) 



which is slightly more complicated than the TBC mixing matrix [14] . but in a better agree- 
ment with the data. Hence it can be considered as an improved TBC mixing. From Eq. (H 
a simple relation between sin 823, sin #13 and the CP phase ip can be derived: 



sin 2 9 



23 



1 / 1 x A 3 

-(1 + A cosy? H cosy?). 

2 8 



where the A 3 term can be neglected. The Jarlskog invariant is given by 

J~ -(Asiny?)/6 + £>(A 3 ). 

Note that ip is not the CP phase in the standard parametrization, but their difference is 
very small, as will be shown in the next section. 

In the case with vanishing CP phase, substituting \J 2 sin 2 #13 for A and using the global-fit 
value for sin 2 # 13 given in Eq. ([TJ , one has 



sin 2 fl 



12 



0.295+°;°°*, sin 2 # 23 ~ 0.607 +a °° 5 



-0.006 
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that agree with the data within la range. Note that by switching A — > —A or ip — > ip + tt, 
sin 2 623 can be brought into the la range in the first octant given in Eq. (CQ). For non- 
vanishing phase, one has 

0.39 < sin 2 #23 < 0.61 
which is within the 3a experimental range 121 ] . 



In addition, it is interesting to see that it may also be regarded as a variant of the 
TM2 trimaximal mixing [ijj] with the TM2 condition \{Uv)a2^ = 1/3 (a = e, fi, r) being 
perturbed by small corrections at the order of 0(A 2 ). One may substitute Ac for A, which 
leads to 

sin 2 9i2 ~ 0.29, sin 2 6 23 ~ 0.61, sin 2 13 ~ 0.025 

which also agree with the data within la range. Hence we refer the mixing matrix given by 
Eq. (|T4|) as trimaximal-Cabibbo mixing matrix. 

One may prefer to perturb £7tbm from the left, i.e. U u = TUtbm- As discussed in Sec. 2, 
in this case, one can move the perturbation matrix to the right by taking transpose. Since 
the discussion is similar, we just give the result. From Eqs. ([6]) we find that the perturbation 
parameters x and y should satisfy 



\—x + y I ~ 



A^ 



t M 2 + M 2 

1 X 



One can verify that 



2 



5A 2 



x 



y 



~ 1 - 3A' 



Xe-* 5A 2 
+ — 



2 8 2 
satisfy the above equations up to higher order terms. Then from Eqs. 



one has 



U ~ 



/ 



V 



J_ 1- 3Ai 
I 1 2 



1 (i_ Ae *_a£) ^{l+ x ^ 



V6 



11A 2 



2 



which leads to a nearly maximal #23, i.e. 

sin 2 6 1 — ~~ 



, 3Ai 
" r 4 



+ 



.J_ _ 3Ai 



1 



1 - 



0(\ 3 ). (15) 



/ 



1 



7 23 



(1 



X 2 

. -). 

Although acceptable, it does not fit the data as well as the trimaximal-Cabibbo mixing 
derived above. One may improve that by adjusting x and y, which is possible for this case, 
but the choice of perturbation parameters is not as straightforward as the one made above, 
so we will leave that for future considerations when more experimental data are available. 
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V. SUMMARY AND DISCUSSIONS 



In this paper a method to construct perturbed mixing matrices and its applications to 
tri-bimaximal mixing are discussed. We show that mixing matrices similar to tri-bimaximal- 
Cabibbo mixing can be constructed by straightforward choices of perturbation parameters. 
However, these mixing patterns fit the data only marginally. To improve that without 
increasing the number of free parameters, we introduce a phenomenological relation, i.e. 
I (^)ei I 2 — I (^)e3 I 2 = 2/3. Two mixing matrices satisfying this relation are constructed. 
The one referred as trimaximal-Cabibbo mixing provides a good two-parameter description 
of the present mixing data and hence can serve as a useful parametrization as long as future 
experimental data does not move the current global fit significantly. Below we discuss briefly 
the phenomenological applications of trimaximal-Cabibbo mixing and its possible theoretical 
explanations. 

As an interesting application, we consider the neutrino mixing probabilities for phase- 
averaged propagation with oscillation phase (Am 2 )L/4:E ^> 1. The explicit formalism has 



already been derived in 



1* 



for triminimal parametrization in which the deviations from 



TBM is described by a perturbation matrix between the (2, 3) and (1, 2) rotation matrices 



that constitute the TBM matrix. By using the formalism given in 18[, the triminimal 
parameters for trimaximal-Cabibbo mixing are found to be 

5A 2 A A 

One also needs the Dirac CP phase in the standard parametrization: 

cosifc ~ cos if — -A 2 cos if sin 2 ip. 

For phenomenological applications, the second term can be ignored since it is at the order of 
O(Xq) at most and hence one has f£, = if. The expressions for flavor mixing probabilities 
can then be derived: 

P VM = 5/9, P v ^ = (7 + 2e + 3e 2 ) /18, P Vr „ Vr = (7 - 2e + 3e 2 ) /18 
P w = (2 - e) /9, P V ^ VT = (2 + e) /9, P w = (7 - 3e 2 ) /18 

where e = Xcosf. The ratio ^^/^ T of the flux to the v T flux arriving at earth which 
measures the violation of /i — r symmetry can be written as 

<V$ T = 1 + 26e 2 /9 
11 



As we can see, the expressions for mixing probabilities are considerably simplified. For more 
details about mixing probabilities, see [18 1 and references therein. 

As to its theoretical background, we note that there are many ways to introduce sub- 
leading contributions. For example, deviations from a LO mixing matrix can be induced 
by radiative corrections, charged lepton corrections, contributions from addition fields and 
terms, etc. In particular, as shown in Sec. 4, the trimaximal-Cabibbo mixing matrix is a 
TM2 trimaximal mixing matrix perturbed by a relatively small correction at the order of 
C(A^-), so in principal it can be explained by trimaximal-mixing models (see, e.g. 
plus an additional small correction. For instance, it is shown in |27[ that a A4 model with 
a 1' (and/or a 1") flavon can lead to the effective neutrino mass matrix: 



M,° = U 



TBM 



a + c — 



V 



fd 





a + 3b + c + d 

a 



^d 




\ 



V 



TBM 



(16) 



c+ij 



where a, b, c and d depend on model parameters. Obviously it leads to a TM2 trimaximal 
mixing. M° can be decomposed as 



^ 1 (A 
1 

v ooi, 



+ 6 



A 1 A 

1 1 1 

V 1 1 V 



+ c 



1 
1 

V 010 / 



\ / 

+ d 



A 

1 



(17) 



It is suggestive to compare it with the one obtained from trimaximal-Cabibbo mixing, i.e. 



mi + 4-ZA21A — ZA 31 A 



M, 



TMC 



TBM 



\ 



V 



m 2 








m 3 



TBM 



(18) 



where rrii are neutrino masses and Ajj 



2 ">sj." w ""j 4—31' 

rrii — rrij. One finds that a (1, 2) off diagonal term 
in the TBM basis is needed. A vanishing CP phase is assumed for simplicity. Note that it 
can also be decomposed as 

x o o x 

1 



M,™ c = M9. + e 



(19) 



V 010 / 
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which can be written as 



M 



TMC 



TBM 



2 ~ t ~ 3 



V2, 



-^e a + 36 + c + d + f 



V 



a - c + | - e 



TBM 



(20) 



where the parameters involved depend on neutrino masses and A. Note that the second 
term in Eq. ( II 9p can be replaced by a combination of ^ o ? 8 j , y 1 1 l j , ^ooij and ^ o 8 S V 

/0 1 l\ (\ 0\ 

e.g. loo or ooo. 
Viooy Voooy 

Such an additional term can be introduced by invoking higher order corrections that 
contribute dominantly to, e.g. the (23) element of the mass matrix, which may not respect 
the Aj symmetry. Alternatively, it can also be achieved by, e.g. adding Higgs triplets (see, 



e.g. 



281]) . Nevertheless, to suppress unwanted contribution to the (23) off diagonal term in 



this tentative model, certain amount of fine-tuning is expected. It would be interesting to 
see a more concrete model that can lead to trimaximal-Cabibbo mixing naturally. We will 
leave that for future work. 
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